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Abstract 

A  computer  model  is  presented  which  describes  the  desorption  of 
liquid  chemicals  of  intermediate  volatility  from  vehicle  interiors.  The 
initial  contamination  level  may  be  derived  by  measurement  of  the  con¬ 
centration  of  the  desorbing  chemical  vapor  as  a  function  of  time.  The 
model  has  been  applied  to  both  allcyd  and  polyurethane  painted  vehicles 
and  used  to  determine  the  residual  contamination  remaining  in  vehicles 
after  decontamination  with  hot  air. 

Introduction 

Chemical  spills  are  not  uncommon  in  modem  society  where  large  volumes  of 
industrial  chemicals  are  transported  by  rail  or  by  road.  Considerable  effort 
is  expended  to  clean  up  spills  from  accidents  and  to  decontaminate  trans¬ 
port  vehicles,  individuals  and  the  surrounding  terrain.  We  have  previously 
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described  studies  detailing  the  removal  of  chemical  contamination  from  ve¬ 
hicle  exteriors  by  natural  weathering  with  sun  and  wind  [1].  However  it  is 
unlikely  that  such  processes  would  be  effective  in  the  enclosed  interiors  of 
vehicles  into  which  contamination  has  been  either  spilled  or  transferred  from 
contaminated  people. 

Personnel  carrying  out  decontamination  frequently  must  wear  protective 
clothing  and  respirators  or  self  contained  breathing  apparatus  as  protection 
against  toxic  vapors.  This  protective  equipment  imposes  a  severe  physiolog¬ 
ical  burden  on  workers  due  to  heat  stress  and  limits  the  effective  operating 
time  for  an  individual.  The  problem  is  exacerbated  in  warm  climates  where 
the  time  which  can  be  spent  in  protective  clothing  while  undertaking  physical 
work  is  greatly  reduced.  Any  process  which  helps  in  the  removal  and  reduc¬ 
tion  of  contamination  will  reduce  the  requirement  for  decontamination  and 
also  will  reduce  the  time  during  which  workers  must  wear  protective  cloth¬ 
ing.  Removal  of  contamination  from  vehicle  interiors  will  reduce  the  hazard 
to  emergency  personnel.  It  is  difficult  to  devise  a  treatment  process  for  chem¬ 
ical  neutralization  that  will  destroy  chemicals  absorbed  into  organic  materi¬ 
als.  However,  physical  processes  are  well  suited  for  the  decontamination  of 
vehicles  and  equipment  that  contain  a  wide  variety  of  materials  ranging  from 
metals  to  glass  and  paints  to  polymers.  The  former  may  be  highly  resistant  to 
chemical  penetration  while  many  of  the  latter  are  highly  absorbent. 

Recent  emphasis  in  removing  chemicals  has  been  on  physical  processes 
that  do  not  involve  aggressive  chemicals  and  are  relatively  benign  to  equip¬ 
ment.  In  particular,  procedures  involving  steam  and  hot  air  have  been  iden¬ 
tified  (2]  as  of  potential  utility.  Higher  temperatures  will  accelerate  removal 
of  contamination  from  surfaces  by  evaporation  and  from  the  bulk  material 
through  desorption.  The  only  limitation  is  the  thermal  stability  of  the  ma¬ 
terials.  Altmann  [3]  has  shown  that  the  temperature  of  hot  air  should  not 
exceed  ISO  °C  if  deterioration  of  thermoplastics  is  to  be  avoided.  For  exam¬ 
ple,  decontamination  processes  for  natural  and  butyl  rubbers  may  be  limited 
to  130-150  °C,  the  temperature  above  which  their  thermal  stability  is  com¬ 
promised. 

The  diffusion  of  a  chemical  from  a  surface  has  been  demonstrated  to 
reveal  a  number  of  interesting  characteristics.  The  kinetics  are  not  always 
strictly  linear;  indeed,  the  level  of  the  chemical  seems  to  play  a  role  in  defin¬ 
ing  its  release  at  low  levels  and  for  particular  simulant  types.  The  release 
mechanism  seems  to  have  a  dependency  on  the  surface  type  from  which  the 
chemical  is  desorbing.  These  facts  coupled  with  the  experimental  complexi- 
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ties  associated  with  monitoring  the  release  and  subsequent  diffusion  of  chem¬ 
ical  vapor  must  be  considered  if  a  single  model  is  to  be  developed  to  describe 
observations. 

We  have  recently  described  such  a  model  that  describes  the  desorption  of 
a  liquid  chemical  from  vehicle  exteriors  [1].  This  model  has  enjoyed  consid¬ 
erable  success  in  its  application  to  accounting  for  simulant  levels  observed  as 
a  function  of  time  and  released  from  a  variety  of  surface  types.  In  the  present 
study  we  show  that  the  model  also  can  be  used  to  describe  the  release  of  chem¬ 
icals  of  intermediate  volatility  from  vehicle  interiors.  As  a  consequence,  the 
initial  contamination  level  may  be  derived  by  measurement  of  the  concentra¬ 
tion  of  the  desorbing  chemical  vapor  as  a  function  of  time. 

The  model  is  of  general  application  to  both  alkyd  and  polyurethane- 
coated  vehicles.  We  have  used  it  to  determine  the  residual  contamination 
of  two  chemical  simulants,  methyl  salicylate  (MS)  and  dimethyl  methylphos- 
phonate  (DMMP),  remaining  in  vehicle  interiors  after  their  decontamination 
with  hot  air.  The  efficacy  of  hot  air  as  a  practical  process  for  removing  chem¬ 
ical  contamination  is  also  demonstrated. 

Modeling  Studies 

There  are  essentially  two  distinct  approaches  to  the  development  of  a  model 
to  describe  the  response  of  a  system: 

1.  a  mechanistic  approach  and 

2.  an  empirical  approach. 

If  the  basic  science  underlying  the  system  whose  response  is  to  be  described 
is  both  well  understood  and  not  too  complex,  a  mechanistic  approach  tends 
to  be  favored.  It  will  yield  a  model  with  gross  properties  as  reflected  in  the 
data  and  may  admit  the  analysis  of  conjectures  in  regard  to  the  system  being 
studied. 

On  the  other  hand,  complex  systems  for  which  accurate  calculations  may 
be  performed  may  be  best  developed  using  an  empirical  approach.  The  dis¬ 
tinction  between  the  two  approaches  is  that  one  is  driven  by  an  exploratory 
motivation  (the  mechanistic  approach)  and  the  other  by  a  data  fitting  ap¬ 
proach  (the  empirical  approach).  That  is,  one  characterizes  mechanisms,  and 
the  other  characterizes  data. 

We  have  developed  an  empirical  approach  to  the  analysis  of  simulant  des¬ 
orption  data  because  of  the  intrinsic  complexity  of  chemical  desorption  pro- 
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cesses  resulting  from  the  wide  variety  of  surfaces  and  materials  present  in 
vehicle  interiors.  While  it  is  possible  to  develop  a  mechanistic  model  of  the 
desorption  of  a  simulant  from  a  simple  flat  painted  surface,  such  an  approach 
is  impractical  for  complex  systems  such  as  vehicles  that  include  a  wide  variety 
of  materials  and  surfaces. 

Previous  studies  have  shown  that  the  evaporation  of  MS  from  alkyd 
painted  plates  followed  first-order  kinetics  and  that  approximately  20%  of 
the  simulant  is  absorbed  strongly  by  the  surface.  The  absorbed  simulant  will 
desorb  more  slowly  creating  a  prolonged  vapor  hazard.  However,  the  evapo¬ 
ration  of  MS  from  polyurethane  paint  could  not  be  described  adequately  by 
the  first-order  model;  the  data  reflected  predominantly  a  zero-order  kinetic 
evaporative  loss  profile.  While  painted  surfaces  form  a  large  part  of  the  inte¬ 
riors  of  Army  Landrovers,  they  are  by  no  means  the  dominant  surface.  Large 
areas  of  vinyl  and  canvas  also  exist  which  will  have  different  rates  and  possibly 
different  mechanisms  of  simulant  evaporation  and  desorption. 

Our  preliminaty  efforts  to  model  mechanistically  the  desorption  data  [1] 
from  vehicles  were  unsuccessful  due  to  the  complexity  of  the  systems.  How¬ 
ever  the  empirical  analysis  approach  proved  to  be  most  valuable.  Although 
an  empirical  model  must  be  based  upon  data,  it  will  be,  if  it  is  to  be  success¬ 
ful,  based  upon  both  the  gross  features  of  the  data  and  the  micro  features  of 
the  data.  The  gross  features  of  the  data  arc  the  qualitative  features  that  we 
notice  recurring  as  we  glance  through  plots  of  replicated  responses  or  glean 
from  discussing  the  observations  with  experimenters.  The  micro  features  of 
data  reflect  variations  from  experiment  to  experiment  and  may  encapsulate 
data  aspects,  which,  if  not  absorbed  into  the  model,  may  compromise  the  ca¬ 
pacity  of  the  model  to  describe  the  data’s  gross  features. 

Gross  features  of  the  data  include  such  traits  as  turning  points  (maxima 
and  minima),  asymptotes,  flatness  of  peaks,  behavior  at  short  times,  inflex¬ 
ions  and  linearizing  transformations.  Micro  features  of  the  data  are  largely 
evident  through  their  interference  with  the  consistency  of  the  gross  features 
and  in  areas  in  the  data  susceptible  to  system  error. 

Our  investigation  of  desorption  studies  has  revealed  the  following  gross 
features: 

1.  an  approximately  exponential  response  phase  followed  by  asymptotic 
leveling, 

2.  a  consistent  small  deviation  from  a  simple  exponential  response,  and 

3.  a  micro  feature  that  implies  that  early  observations  may  be  susceptible 
to  data  collection  error. 
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Assembling  these  into  a  simple  empirical  model  has  yielded  the  following 
function  which  we  have  found  to  be  able  to: 

1.  describe  consistently  the  simulant  desorption  data,  and 

2.  highlight  experimental  problems  (e.g.,  contamination)  which  will  thus 
draw  attention  to  data  sets  that  are  not  fully  consistent  with  the  model 
proposed. 

y  =  P,(  (1) 

where  y  is  the  expected  cumulative  dosage.  For  this  function  we  see  that 

Pi  is  the  asymptotic  level  of  desorption  data, 

P2  is  an  estimate  of  the  desorption  rate  constant, 

P3  is  a  measure  of  the  degree  of  adherence  to  linear  kinetics,  and 
P4  is  an  offset  parameter  which  serves  to  accommodate  uncertainties  in 
the  actual  sample  time  particularly  associated  with  early  observations. 

Model  Parameters 

If  used  systematically  in  a  data  fitting  sense  eq  1  has  been  found  to  yield 
constant  and  sensible  estimates  for  its  parameters.  The  model  describes  the 
desorption  of  moderate  levels  of  contaminant  which  do  not  reach  saturated 
vapor  levels.  It  will  not  apply  for  very  small  or  very  large  values  of  surface 
contamination.  Furthermore,  if  we  require,  calculates  and  their  errors  can 
be  accurately  determined.  In  this  section  of  the  report  we  outline  a  standard 
approach  to  model  fitting. 

The  Asymptotic  Value:  P, 

A  major  aim  of  the  investigation  has  been  to  estimate  the  model  parameter 
Pi,  which  is  experimentally  related  to  the  asymptotic  value  of  the  cumulative 
dosage  of  vapor  released  by  the  vehicle,  and  to  calibrate  Pi  as  an  indicator 
of  initial  exposure  to  simulant.  The  latter  addresses  the  question  of  “given 
a  profile  of  simulant  desorption  through  time  can  we,  by  fitting  eq  1  to  this 
profile,  infer  the  initial  exposure  of  the  surface  from  which  the  simulant  des¬ 
orbed”.  This  issue  was  addressed  in  an  earlier  report  [1].  Flattening  of  the 
response  will  be  evident  from  the  horizontal  tendency  of  the  data  when  plot¬ 
ted  in  either  linear  or  semi-logarithmic  mode.  Indications  of  problems  in  this 
regard  are  evident  from  a  final  estimate  of  Pi  that  is  much  larger  than  the 
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final  datum  value,  an  inaccurate  estimate  (error  in  excess  of  20%)  and  slow 
convergence  to  the  final  fit. 

The  Time  Offset:  P4 

Although  only  a  micro  feature  of  the  data,  the  capacity  of  P, « to  (a)  improve 
the  fit  and  (b)  compromise  the  fitting  procedure  is  quite  dramatic.  The  time 
offset  dominates  the  first  unit  (1  h  or  so)  of  eq  1,  though  if  it  is  not  carefully 
set  or  controlled  it  may  drift  to  a  value  larger  than  the  time  point  coinciding 
with  the  first  observation  and  cause  a  numerical  error. 

The  Desorption  Rate  Constant:  Pa 

For  the  moment,  neglecting  P»  and  setting  P3  equal  to  1,  we  see  that  Pi  defines 
the  rate  of  rise  of  eq  1. 

Writing  y  *  ln(l  -  y/Pi)  (2) 

we  find  that  y  «  -Pit  +  log(l/Pi)  (3) 

Thus  we  see,  using  observations  after  three  units  or  so,  that  initial  es¬ 
timates  for  Pi  can  be  obtained  from  the  slope  of  the  semi-logarithmic  plot 
given  by  eqs  2  and  3.  The  rate  of  rise  of  eq  1  is  defined  by  P2.  Its  value  seems 
ordinarily  to  adjust  to  around  0.1/h  for  MS  and  0.2/h  for  DMMP.  The  most 
sensitive  domain  of  y  to  P2  is  from  1  h  to  3-4  h.  Whilst  improved  estimates 
of  Pi  will  be  achieved  using  the  inverse  observation  form  of  data  weighting, 
this  is  not  recommended  as  we  are  usually  motivated  primarily  to  explain  the 
higher  observations  (i.e.,  asymptotic  values)  and  secondarily  to  explain  the 
overall  response.  In  any  case,  there  are  rarely  any  problems  with  the  estima¬ 
tion  of  Pi. 

The  Deviation  from  Linear  Desorption:  P3 

Our  earlier  investigations  have  shown  that  under  certain  conditions  and  at 
certain  stages  during  the  desorption  process,  the  kinetics  of  desorption  be¬ 
come  nonlinear.  The  purpose  of  P3  is  to  facilitate  that  trait  in  the  model. 
Differentiating  and  rearranging  eq  1,  we  find  that 

y  -  -P2j/(f,P,"P3(l  -  P\fy)]  (4) 

where  t'  =  f  -  P4 

which  reduces  to  a  first  order  linear  differential  equation  if  P3  =  1.  If  P3 
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encurses  to  a  boundary  then  usually  we  find  that  contamination  or  some  other 
problem  exists  in  the  data.  Otherwise,  based  on  response  measurements  in 
its  dominant  range  (4-1  Sh),  Pj  is  usually  consistently  and  sensibly  estimated. 


Sensitivity  Analysis 

In  our  discussion  above  on  the  model  parameters,  we  alluded  to  the  “domi¬ 
nant  ranges”.  In  particular  we  suggested  that  the  most  sensitive,  or  dominant 
ranges  of  the  response  function  eq  1  to  P2  and  P3  were  1  h  to  3-4  h  and  4  h 
to  15  h,  respectively.  For  any  response  function  y  with  parameter  6, 


dy‘  = 


r  =  {o\,o\  -  K) 

.  djLdB.  .  .  • . 

d0\  d°x  +  d0\  d°2  d6'nde" 


dy ’  =  Sdd 


•t 


(5) 


where 

dr  =(doid<r2,.-,d0-n) 

and 

*  (d£  d£_  d£\ 

\d9\‘  d0'2'  "  '  dd'J 

is  the  dynamic  sensitivity  of  y  to  0  which  can  be  calculated  by  post-multiplying 
both  sides  of  eq  5  by  the  identity  matrix  L  with  dimension  n  x  n.  Then  we  have 


dy'L  =  SdO"L 
dy*  =  SdO* 

and  post-multiplying  by  the  inverse  gives 

5  =  dy*(d0*)~l 

where  each  element  in  S  is  given  by: 

O  _  dy*  _  dy/dyo  _  dlny^ 
'  d0;  ~  d0i/9iO  d  In  9\ 


SENSITIVITY 
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TIME  (h) 

Figure  1:  Sensitivity  range  of  model  parameters. 


The  dynamic  sensitivity  is  of  considerable  importance  for  two  reasons: 

1.  it  helps  to  isolate  parameters  in  y  which  may  be  difficult  to  estimate  and 

2.  it  helps  to  ensure  that  the  distribution  of  observations  with  time  is  ap¬ 
propriate  to  allow  accurate  estimation  of  all  the  model's  (response)  pa¬ 
rameters. 

In  Figure  1  we  present  a  plot  of  the  sensitivity  5  for  each  of  the  parameters 
in  eq  1  and  the  following  observations  can  be  made: 

1.  The  dominant  range  of  Pa  is  that  associated  with  the  earliest  values  of 
y  (i-c.,  t  <  1  h). 

2.  The  dominant  range  of  Pi  is  from  1  h  to  3-4  h. 

3.  The  dominant  range  of  Pi  is  from  4  h  to  15  h. 

4.  The  dominant  range  of  P\  is  for  time  values  beyond  IS  h. 

5.  The  early  time  range  is  dominantly  influenced  by  two  parameters  Pi, 
and  Pa,  however  because  Pi  also  dominates  a  later  response  region  its 
resolution  in  association  with  the  observations  there  will  permit  the  res¬ 
olution  of  Pa  in  the  early  response  range. 
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6.  Because  oi  the  unique  time  ranges  in  which  each  of  the  parameters  has  a 
major  effect  on  the  response,  accepting  (5),  an  adequate  sampling  strat¬ 
egy  out  to  48  h  should  characterize  the  response,  and  hence  provide  an 
adequate  set  of  observations  for  the  identification  of  each  parameter 
Our  current  sampling  strategy  includes  observations  at  the  following 
time  points:  1(1)8,  10(2)24,  27(3)36, 40(4)48,  where  m(n)p  means  ad¬ 
vancing  from  m  h  to  p  h  in  steps  of  n  h. 

Experimental 

Desorption  Chamber 

The  chamber  in  which  the  studies  were  carried  out  was  constructed  of  stain¬ 
less  steel  and  was  thermally  insulated  and  temperature  controlled.  Tbtal  vol¬ 
ume  of  the  chamber  was  93  m3.  Mechanical  circulation  of  the  air  inside  the 
chamber  ensured  rapid  mixing  of  desorbing  vapor. 

Our  toxic  chemical  simulants,  DMMP  and  MS,  were  chosen  because  of 
their  intermediate  volatilities.  They  are  commonly  used  as  simulants  for  other 
toxic  industrial  and  military  chemicals  for  this  reason. 

Analysis 

The  MS  vapor  concentrations  were  monitored  by  sampling  into  propylene 
glycol  in  sequential  sampler  bubblers.  The  sequential  samplers  each  held 
12  bubblers  containing  5  mL  propylene  glycol  as  the  absorbing  medium  and 
were  programmed  to  sample  the  chamber  over  a  48-h  time  span.  The  bubbler 
contents  were  analyzed  subsequently  for  MS  by  UV  spectroscopy.  Similarly, 
DMMP  was  sampled  into  diethyl  phthalate  in  sequential  sampler  bubblers. 
Bubbler  contents  were  analyzed  for  DMMP  by  gas  chromatography  using  a 
flame  photometric  detector. 

Vehicles 

The  vehicles  chosen  were  Army  Landrovers  complete  with  canvas  canopies. 
The  vehicles  were  coated  with  two  paint  systems;  in-service  Matt  Olive  Drab 
Alkyd  paint  or  Low  Gloss  Olive  Drab  Aliphatic  Polyurethane  paint  (PUP)  to 
GPC-P  154/3  [4].  The  vehicle  interiors  essentially  consisted  of  painted  metal 
on  the  floors  and  dashboard,  canvas  in  the  canopy  and  vinyl  seat  covers. 


2310 


AMOS,  LEAKE,  AND  BOSTON 


For  calibration  purposes,  known  quantities  of  liquid  simulant  were  ap¬ 
plied  as  5-/iL  droplets  from  a  syringe  to  the  various  interior  surfaces  of  the  ve¬ 
hicle.  For  the  hot  air  decontamination  studies,  liquid  simulants  were  sprayed 
on  to  the  interior  surfaces  in  known  quantities  from  a  syringe. 

The  contaminated  vehicle  was  driven  into  the  sealed  chamber  which  was 
temperature  controlled  at  25  ±  1  °C.  Simulant  desorbing  inside  the  vehicle 
was  monitored  by  collection  for  later  analysis  in  four  sets  of  sequential  sam¬ 
plers  arranged  at  head  height  in  the  driving  and  passenger  positions.  In  the 
studies  on  decontamination  with  hot  air,  the  contaminated  vehicles  were  al¬ 
lowed  to  stand  for  30  minutes  to  permit  simulant  to  soak  into  materials  and 
to  mimic  more  closely  a  real  situation.  They  were  then  decontaminated  with 
hot  air  at  130  °C  for  either  30  or  60  minutes.  Appropriate  hot  air  at  a  flow 
rate  of  50  L/min  was  supplied  from  a  Karcher  FB  60  E  Hot  Air  Generator 
and  blown  into  the  rear  compartment  through  two  flexible  ducts.  The  ve¬ 
hicles  were  placed  in  the  chamber  immediately  once  decontamination  was 
complete.  Analysis  of  all  data  about  the  desorption  of  simulants  from  inside 
the  vehicles  was  carried  out  using  CONSAM  which  is  the  interactive  version 
of  the  SAAM  modeling  program  and  allows  the  user  to  develop  mathematical 
models  to  flt  experimental  data  [4]. 


Results  and  Discussion 

Previous  studies  have  shown  that  the  evaporation  of  MS  from  aikyd  painted 
surfaces  follows  flrst-order  kinetics  and  that  the  evaporation  of  MS  from 
polyurethane  paint  shows  a  zero-order  kinetic  loss  profile.  That  simple  first- 
order  kinetics  do  not  describe  the  evaporation  of  simulants  from  both  aikyd 
and  polyurethane  vehicles  is  evident  from  considerations  of  the  nature  of 
change  of  the  surface  area  of  the  simulant  as  it  evaporates.  It  is  likely  that 
the  nature  of  simulant  kinetics  will  be  initially  zero  order.  Then,  as  the  sim¬ 
ulant  evaporates  to  a  stage  where  its  surface  area  is  reducing  linearly,  the 
kinetics  will  be  first-order,  and  finally  nonlinear.  Other  absorbent  materials 
present  inside  vehicles,  such  as  canvas  and  vinyl,  may  complicate  the  kinetics 
of  simulant  desorption.  Chemical  contamination  can  remain,  held  by  cap¬ 
illary  action  in  cracks  and  crevices  even  when  the  major  part  of  the  vehicle 
surface  is  contaminant-free.  Regardless  of  the  nature  of  the  surface,  outward 
diffusion  of  contamination  vapor  is  accelerated  by  the  increased  temperatures 
resulting  from  decontamination. 
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Figure  2:  Desorption  of  DMMP  from  [A]  PUP  and  (B]  ALKYD  vehicle  inte¬ 
riors  -  experimental  data  and  CONSAM  fitted  equations.  The  unbroken  line 
is  the  CONSAM  fit. 


There  is  no  way  of  determining  the  various  phases  of  desorption  and  it 
was  more  useful  to  explore  an  empirical  rather  than  a  mechanistic  approach 
as  the  desorption  profile  may  span  several  regions  of  applicability  of  a  mecha¬ 
nism.  Furthermore,  we  are  more  concerned  with  the  accuracy  of  a  prediction 
and  the  estimation  of  the  residual  contamination  than  with  desorption  mech¬ 
anisms.  We  therefore  attempted  some  empirical  experimentation  with  the 
model  developed  for  vehicle  exteriors  which  was  potentially  flexible  enough 
to  absorb  the  key  kinetic  features  alluded  to  above  and  to  explain  the  data  for 
both  types  of  vehicles. 


Modeling  of  Simulant  Desorption  from  Vehicles  Interiors 

Evaporation  of  MS  simulant  from  the  painted  external  surfaces  of  the  ve¬ 
hicles  has  been  shown  previously  to  be  first-order  for  ALKYD  vehicles  and 
zero-order  for  PUP-painted  vehicles.  The  interiors  of  both  vehicles  have  con¬ 
siderable  surface  areas  of  canvas  and  vinyl  which  may  desorb  with  first-order 
kinetics  and  tend  to  mask  the  behavior  of  the  polyurethane  paint.  Further¬ 
more,  the  constitution  and  gloss  level  of  the  painted  surface  has  a  consider- 
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Figure  3:  Desorption  of  MS  from  [A]  PUP  and  [B]  ALKYD  vehicle  interiors 
-  experimental  data  and  CONSAM  fitted  equations.  The  unbroken  line  is  the 
CONSAM  fit. 


able  influence  on  the  desorption  and  evaporation  kinetics.  However  this  may 
well  be  masked  by  the  large  proportion  of  other  absorbent  materials. 

Simulant  was  applied  to  the  vehicle  interior  surfaces  as  described  in  the 
Experimental  section  and  the  cumulative  desorption  of  simulant  into  the  ve¬ 
hicle  was  measured.  Figures  2  and  3  summarize  the  fits  of  the  above  model  to 
a  series  of  studies  relating  to  three  levels  of  DMMP  and  MS  simulant  contam¬ 
ination  density,  respectively,  for  each  of  the  two  types  of  vehicles.  Each  plot  is 
the  mean  of  four  replicates.  The  vehicles  were  not  subject  to  decontamination 
in  this  series  and  vapor  sampling  started  immediately  after  contamination. 

The  fitted  plots  were  produced  from  the  experimental  data  using  CON¬ 
SAM  and  the  developed  empirical  model,  eq  1.  The  data  were  used  to  define 
the  model  and  to  relate  initial  contamination  levels  to  parameters  derived 
from  the  model.  The  derived  parameters  and  associated  standard  deviations 
are  given  in  Thble  1  for  DMMP  and  in  Thble  2  for  MS. 

Of  particular  interest  is  the  fact  that  Pj  «  1  for  MS  on  both  vehicles, 
implying  that  MS  desorption  is  linear.  In  contrast,  Pi  «  0.5  for  DMMP, 
implying  that  desorption  departs  considerably  from  linearity. 

With  the  initial  conditions  (initial  simulant  contamination  level)  known 
for  each  experiment,  the  estimated  P\  for  studies  in  which  neither  Pi  nor 
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Tfcble  1:  Derived  model  parameters  for  DMMP  desorption  from  vehicle 


interiors. 


Veh.a 

DM6 

Pi 

P2 

Pi 

P4 

Pe 

858 

28.7  ±  4.0 

0.179  ±0.022 

0.430  ±  0.026 

0.894  ±  0.004 

Pe 

481 

11.0  ±0.3 

0.220  ±  0.005 

0.586  ±  0.018 

0.844  ±  0.040 

Pe 

309 

7.2  ±  0.6 

0.296  ±  0.021 

0.383  ±  0.022 

0.932  ±  0.020 

Ad 

858 

27.0  ±4.1 

0.118  ±0.016 

0.452  ±  0.018 

0.895  ±  0.028 

Ad 

550 

13.0  ±  0.6 

0.235  ±  0.008 

0.490  ±  0.018 

0.954  ±  0.017 

Ad 

309 

8.0  ±  0.6 

0.154  ±0.007 

0.591  ±  0.037 

0.795  ±  0.099 

“Vehicle  type. 

bDimethyl  mcthylphosphonaic  in  milligrams.  ePUP. 

dALKYD. 

Ihble  2: 
riors. 

Derived  model  parameters  for  MS  desorption  from  vehicle  intc- 

Veh.a 

MS* 

Pi 

Pz 

Pi 

P. 

Pe 

880 

110  ±0.1 

0.106  ±  0.003 

0.902  ±  0.015 

0.558  ±  0.006 

?e 

598 

5.6  ±  0.04 

0.114  ±  0.007 

0.958  ±  0.003 

0.851  ±  0.102 

?e 

422 

4.1  ±  0.02 

0.106  ±0.007 

1.055  ±  0.028 

0.591  ±0.105 

Ad 

880 

11.1  ±0.11 

0.137  ±0.006 

0.817  ±0.022 

0.633  ±  0.084 

Ad 

598 

7.9  ±0.03 

0.113  ±0.005 

0.961  ±  0.017 

0.692  ±  0.067 

Ad 

422 

3.5  ±  0.01 

0.084  ±  0.007 

1.150  ±0.030 

0.410  ±0.119 

“Vehicle  type.  kMethyl  salicylate  in  milligrams.  CPUP.  dALKYD. 


P3  was  constrained  was  regressed  on  simulant  dosage,  depicted  in  Figure  4. 
Within  the  region  under  study,  P\  is  related  to  the  initial  contamination  level 
as  follows: 

Pi  =  38.0(/C)  -  5.0  for  DMMP  (6) 

and  Pi  =  17.0(/C)  -  3.0  for  MS 

where  IC  is  the  initial  contamination  level  in  grams. 

The  desorption  rate  constant,  P2,  has  a  value  of  around  0.1/h  for  MS  and 
0.2/h  for  DMMP.  This  constant  may  be  influenced  by  several  physical  proper¬ 
ties  of  the  simulants.  In  general  terms,  surface  tension  will  affect  the  spread- 
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Figure  4:  Plot  of  Pi  vs  initial  contamination  level  of  simulants  [A]  DMMP 
and  [B]  MS  for  ALKYD  and  PUP  vehicles  using  unconstrained  values  of  P\. 


Thble  3:  Values  of  parameter  P2  and  physical  properties  of  simulants. 


DMMP 

MS 

P2I h 

0.2 

0.1 

Surface  tension  (dyne/cm) 

35.9 

38.8 

Vapor  Pressure  at  25  °C  (torr) 

0.613 

0.165 

Solubility  parameter 

10.6 

10.6 

ing  of  simulant  over  the  surface  and  the  penetration  into  cracks  and  crevices. 
\fepor  pressure  will  affect  the  evaporation  of  the  simulant  from  the  surface. 
The  solubility  of  the  simulants  in  the  various  organic  materials  of  the  vehicles 
will  affect  their  diffusion  into  and  out  of  the  bulk  materials.  Molecular  size 
and  specific  interactions  between  simulant  and  materials  will  influence  the 
rate  of  desorption  from  the  surfaces. 

We  have  detailed  some  relevant  physical  properties  of  the  simulants  in  Ta¬ 
ble  3,  above.  The  lower  surface  tension  and  higher  vapor  pressure  of  DMMP 
compared  to  MS  result  in  increased  spreading  and  evaporation  rates  and  are 
reflected  in  the  values  of  P2.  The  solubility  parameters  for  DMMP  and  MS 
are  the  same,  implying  that  the  simulants  are  equally  absorbed  by  the  organic 
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(a)  6_  W«' 


Figure  5:  Desorption  of  MS  from  vehicle  interiors  (A]  PUP,  [B]  ALKYD  after 
hot  air  decontamination  for  30  or  60  minutes  •  experimental  data  and  CON- 
SAM  Fit  for  various  levels  of  initial  contamination.  The  unbroken  line  is  the 
CONSAM  fit 


materials  in  the  vehicles.  Dimethyl  methylphosphonate  is  a  smaller  molecule 
than  MS  and  will  therefore  diffuse  out  of  absorbent  materials  more  rapidly. 
Specific  interaction  between  MS  and  polymeric  materials,  possibly  due  to  hy¬ 
drogen  bonding,  will  also  reduce  the  relative  rate  of  desorption  of  MS.  Such 
complexity  governing  the  desorption  of  organic  liquids  from  complex  systems 
mitigates  against  the  utility  of  a  mechanistic  approach  to  desorption  and  re- 
infc 'v«s  the  benefits  of  the  empirical  data  analysis  approach. 

Desorption  of  Simulants  from  Contaminated  Vehicles  after 
Decontamination  with  Hot  Air 

Ibtal  cumulative  dosages  of  MS  desorbing  from  PUP  and  ALKYD  vehicle 
interiors  after  decontamination  with  hot  air  for  30  minutes  are  shown  in  Fig¬ 
ure  5.  Similarly  dosages  of  DMMP  desorbing  from  the  vehicle  interiors  after 
decontamination  are  displayed  in  Figure  6.  The  fitted  plots  were  produced 
from  the  developed  empirical  model  using  CONSAM.  Levels  of  residual  sim- 
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Figure  6:  Desorption  of  DMMP  from  vehicle  interiors  (A]  PUP,  [B]  ALKYD 
after  hot  air  decontamination  for  60  minutes  -  experimental  data  and  CON- 
SAM  fit  for  various  levels  of  initial  contamination.  The  unbroken  line  is  the 
CONSAM  fit. 


"fable  4:  Derived  values  of  parameter  P\  and  calculated  values  of  residual 
MS  simulant  after  hot  air  decontamination. 


Vehicle 

Simulant 
Level (g) 

Decontamin. 
Time  (mins) 

P\ 

Residual 
MS  (mg) 

Decontamin. 
Effic.  (%) 

PUP 

1.14 

60 

ND 

ND 

100 

ALKYD 

1.09 

60 

ND 

ND 

100 

PUP 

2.38 

60 

0.95 

250 

90 

ALKYD 

2.24 

60 

5.8 

540 

77 

PUP 

1.10 

30 

3.7 

410 

63 

ALKYD 

1.10 

30 

3.9 

430 

60 

PUP 

2.20 

30 

6.8 

600 

73 

ALKYD 

2.24 

30 

9.0 

730 

67 
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"fable  5:  Derived  values  of  parameter  P\  and  calculated  values  of  residual 
DMMP  simulant  after  hot  air  decontamination. 


Vehicle 

Simulant 
Level  (g) 

Decontamin. 
Time  (mins) 

Pi 

Residual 
DMMP  (mg) 

Decontamin. 
Effic.  (%) 

PUP 

1.19 

60 

3.90 

247 

79 

ALKYD 

1.22 

60 

4.93 

274 

78 

PUP 

2.07 

60 

8.62 

372 

82 

ALKYD 

2.39 

60 

6.52 

316 

87 

PUP 

4.63 

60 

14.4 

526 

89 

ulant,  the  contamination  level  after  decontamination,  were  calculated  from 
eq  6.  "fables  4  and  5  summarize  details  of  the  initial  contamination  levels, 
decontamination  times,  values  of  the  derived  parameter.  Pi,  and  decontami¬ 
nation  efficiency. 

In  these  tables  we  also  present  the  efficiency  of  hot  air  decontamination. 
It  is  evident  that  treatment  with  hot  air  at  130  °C  for  periods  around  60  min¬ 
utes  is  an  effective  way  to  reduce  contamination.  There  is  little  evidence  of 
any  difference  in  performance  between  the  polyurethane  and  alkyd  painted 
vehicles.  This  is  due  to  the  large  area  of  canvas  and  vinyl  which  are  common 
to  the  test  vehicles  and  which  mask  differences  between  the  two  coatings.  Ve¬ 
hicles  in  which  the  interiors  are  mainly  painted  metal  may  decontaminate  at 
different  rates  dependant  on  the  paint  properties.  Methyl  salicylate  is  more 
readily  decontaminated  than  DMMP  even  though  the  volatility  of  the  latter 
is  higher.  This  may  be  attributable  to  the  greater  interaction  and  solubility  of 
DMMP  in  the  paints,  canvas  and  vinyl  seat  material. 

Commercial  transport  vehicles  are  coated  with  a  variety  of  paint  types  and 
have  a  range  of  interiors  fittings  and  materials.  These  interiors  are  likely  to 
be  chemical  absorptive  and  similar  to  the  alkyd  coated  vehicle  in  our  study. 
As  a  consequence,  hot  air  decontamination  may  be  considered  effective  for  a 
wide  range  of  commercial  vehicles.  Decontamination  efficiency  will  depend 
on  the  interaction  of  contamination  with  the  materials  as  well  as  the  volatility 
of  the  chemical  involved.  For  chemicals  with  a  volatility  less  than  or  equal 
to  that  of  MS,  decontamination  with  hot  air  at  130  °C  for  60  minutes  is  an 
effective  process.  The  model  we  have  described  is  flexible  enough  to  include 
such  variations. 
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Notation 

DMMP  dimethyl  methylphosphonate 
MS  methyl  salicylate 

Pi  the  asymptotic  level  of  desorption  data 

P2  an  estimate  of  the  desorption  rate  constant 
Pi  a  measure  of  the  degree  of  adherence  to  linear  kinetics 
P4  an  offset  parameter  which  serves  to  accommodate  uncertainties  in  the 

actual  sample  time  particularly  associated  with  early  observations 

5  dynamic  sensitivity  of  y  to  a  given  parameter 

t  time  from  start 

6  a  generic  parameter 

y  cumulative  dose  of  toxic  chemical 
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